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Possible variations of the quark-meson couphng (QMC) model are examined 

in which the bag constant decreases in the nuclear medium. The reduction 

is supposed to depend on either the mean scalar field or the effective mass 

of the nucleon. It is shown that the electric and magnetic radii of the bound 

;_( ■ nucleon are almost linearly correlated with the bag constant. Using the fact 

C^ , 

that the size of the bound nucleon inside a nucleus is strongly constrained by 

y-scaling data in quasielastic, electron-nucleus scattering, we set a limit for 

the reduction allowed in the bag constant for these two models. The present 

study implies that the bag constant can decrease up to 10-17 % at average 

nuclear density, depending on the details of the model. 
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I. INTRODUCTION 

Whether a nucleon bound in a nuclear medium (precisely speaking, a bound "nucleon" 
is a quasiparticle which has the same quantum numbers as those of the nucleon) would 
alter its properties from those of a free nucleon is of fundamental interest in nuclear physics 
||l|-^. It is a traditional approach that nucleonic (and mesonic) degrees of freedom play an 
essential role in describing nuclei, but it is also true that those nucleons themselves have 
substructure. The idea of employing quark degrees of freedom to describe certain aspects 
of nuclear structure was motivated in part by the discovery of the EMC effect, namely, 
the discovery that the structure functions of the nucleon [as measured in deep inelastic 
scattering (DIS)] depend on the nucleus in which the nucleon is bound ||^. It has been 
widely accepted that DIS is well described by quark and gluon dynamics, i.e., quantum 
chromodynamics (QCD). Thus, it is reasonable to expect that quarks inside a nucleus may 
play important roles in some contexts, and their wave functions could be strongly modified 
by the surrounding nuclear environment. 

There is strong evidence for a reduction of the effective nucleon mass in the nuclear 
medium. In quantum hadrodynamics (QHD) ||l| and the quark-meson coupling model 
(QMC) 10-^, the effective mass of nucleon at normal nuclear density, in mean-field ap- 
proximation (MFA), is respectively predicted to be roughly 40-50% and 20-30% smaller 
than the value in free space. The mechanism for this reduction is the interaction of the nu- 
cleon with the mean scalar field. QCD sum rules independently lead to a similar prediction 
from an apparently different mechanism, i.e., the reduction of the quark condensate as the 
density increases ||^. In contrast, whether the size of the nucleon changes and by how much, 
are relatively less clear. The validity of Brown-Rho scaling for the nucleon radius could 
be much more questionable than it is for hadron masses. 



The successes of QHD leave little doubt that relativistic nuclear phenomenology is es- 
sential in describing the bulk properties of nuclear matter as well as the properties of finite 
nuclei. A compelling feature of relativistic mean-field theory has been the reproduction 
of spin-orbit splittings in finite nuclei and the description of spin observables in nucleon- 
nucleus scattering, after calibrating the model parameters by the equilibrium properties of 
symmetric nuclear matter (at normal nuclear density). These successful results could be 
ascribed to the reduction of the effective mass of the nucleon |jl| . 

To incorporate the composite nature of the nucleon without losing the successful phe- 
nomenology of QHD, Guichon proposed a simple hybrid model (the quark- meson coupling 
(QMC) model) in which nuclear matter is described by non-overlapping, MIT bag nucle- 
ons 0. By analogy with QHD, QMC describes the bulk properties of nuclear systems using 
scalar (a) and vector (u) meson mean fields. However, the nucleon bound in nuclear medium 
here is no longer a point particle, it has substructure - quarks confined inside the nucleon 
bag. It is the quark, rather than the nucleon itself, that is coupled to the a and uj fields 
directly. As a result, the internal structure of the bound nucleon is modified by the sur- 
rounding medium with respect to that in free space. The relative simplicity of the QMC 
model at the hadronic level, together with its phenomenological successes, makes it suitable 
for many applications in nuclear physics [|5i [l0| , [ri| . There are similar models which use the 



harmonic oscillator [|I0| or the color dielectric model []12[ instead of the MIT bag for the 
nucleon. 

Within QMC, the small mass of the quark implies that the lower component of the 
quark wave function will be enhanced rapidly by the change of its environment (as the 
a field strength increases), with a consequent decrease in the scalar baryon density. As 
the scalar baryon density itself is the source of the a field, this provides a mechanism 
for the saturation of nuclear matter where the quark substructure plays a vital role. The 
extra degrees of freedom, corresponding to the internal structure of the nucleon, lead to a 
reasonable value for the nuclear incompressibilty, once the corresponding quark and meson 
coupling constants, g^ and g^, are determined to reproduce the empirical values for the 



saturation density and binding energy of symmetric nuclear matter. One drawback in this 
approach is that the effective scalar and vector fields are significantly reduced, which in turn 
leads to somewhat insufficient spin-orbit interaction in finite nuclei (within MFA) ||13|| . 

With this difficulty in mind, Jin and Jennings postulated that the bag constant (B) might 
decrease along with a reduction of the effective nucleon mass {m*^) |T^. In their model, 
the bag constant had to be reduced dramatically (approximately 35% at normal nuclear 
density) to achieve "the large and cancelling scalar and vector potentials" which arise in 
QHD and lead to the large nuclear spin-orbit forces. One undesirable side effect is that the 
nuclear incompressbility once again tends to be overpredicted (one of the original successes 
of Guichon's model was that it generated a reasonable value for the incompressibility) . In 
addition, the nucleon size was found to increase unpleasantly because of the reduction of 
the bag constant. (For a recent study, see also Ref. [p!5| .) 

On the other hand, there are already strong, experimental constraints on the extent 
of such changes of nucleon size p6|-p!^. Furthermore, new quasi-elastic (e, e') and {e,e'N) 
experiments at facilities such as Mainz, MIT-Bates and especially TJNAF should improve 
these constraints in the near future. 

In this paper we examine this issue, namely, how much the bag constant can be allowed 
to decrease in nuclear matter given those experimental constraints. For this purpose, we first 
solve self-consistently the equations of motion for the quarks and mesons in the modified 
QMC models ||14|. Then we apply the solutions obtained to calculate the electromagnetic 
form factors of the nucleon in medium. From these calculated form factors, we extract 
the root-mean-square (r.m.s.) magnetic- and electric- radii. By comparing these with the 



experimental data for the change of the nucleon size ||16|- |I9|J we can constrain the reduction 

of the bag constant. For simplicity, possible off-shell effects of the bound nucleon in medium 

I9| , |20[] are neglected in the present treatment, but should be added in future, more elaborate 



studies. However, since each nucleon bag contains three valence quarks in the ground state 
which move independently under the spherical cavity approximation to the MIT bag ^ , the 
spurious center-of-mass motion has to be removed. In this study, we use the Peierls-Thouless 



(PT) projection method combined with a Lorentz contraction for the nucleon internal wave 
function, as in Ref. |^. This technique for implementing the center-of-mass and recoil 
corrections has been quite successful in the case of the electromagnetic form factors of 



nucleon in free space [21 



In conventional nuclear physics, the nucleon is immutable and the electromagnetic prop- 
erties of a nucleus are often described by a combination of individual nucleon contributions 
and various meson exchange current (MEC) corrections. Although we believe that the two 
body meson exchange is small in the present case, a full calculation of a nuclear form factor 
requires a consistent treatment of all MEC corrections. Since we will ignore the explicit 
MEC corrections here for simplicity, our results should be viewed as somewhat tentative 
with a need for future elaboration. 

In section II, we briefly summarize our version of the modified QMC model with the 



option for the reduction of the bag constant. As in Ref. [Q, B is related to either a or 
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however, the parameters are quite different because of the improved technical treatment of 
center-of-mass motion. In sec. Ill, we first discuss the experimental status of the constraints 
on the nucleon electromagnetic radii. The basics of the ^/-scaling analysis are presented. 
Then we show the formalism used to calculate the nucleon electromagnetic form factors and 
r.m.s. radii. In sec. IV we discuss pionic corrections and the effects of meson-exchange 
currents. The numerical results are discussed in sec. V and final conclusion is given in 
sec. VI. 

II. ASPECTS OF THE GENERALIZED QMC MODEL 

We first review the standard treatment of the QMC model P,^,p!3| , which does not allow 
for any variation of the bag constant, B, in the nuclear medium. 

In QMC, nuclear matter is described in terms of non-overlapping MIT bags, and the 
quarks inside the bags interact directly with the scalar and vector meson fields, where, the 
motion of the quarks is highly relativistic and typically they move much faster than the 



nucleoli in the medium (the Fermi motion of the nucleon) . Thus it is reasonable to assume 
that the quarks always have sufficient time to adjust their motion as the nucleon responds 



to the medium it is in so that they stay in the lowest energy state ||T3[ . In the instantaneous 
rest frame of a nucleon in symmetric nuclear matter (where the p meson mean-field will 
vanish), the Lagrangian density for the QMC model may be given, 

Cq = q{ij''d^ - mq)q9v - BqOv + c/^gga - glq^i^^quj^ - ^mla^ + -mlu^, (1) 

where mq is the current quark mass, Bq is the bag constant in vacuum (and will remain 
fixed), g'^ and g^ are the corresponding quark and meson coupling constants, and By is a 
step function which is one inside the bag volume and vanishes outside. 

In MFA, the meson fields are treated as classical fields, and the quark field q{x) inside 
the bag satisfies the equation of motion 

[^Y^, - (m, - gla) - g!^uJ^']q{x) = 0, (2) 

where a and uJ denote the constant mean values of the scalar and the time component of 
the vector field in symmetric nuclear matter. The normalized solution for the lowest state 
of the quark is given by [^,|| 



v47r 



jo(a;r/i?) 
^iPqCr ■rji{xr/R) ^ 



0{R-r)x,, (3) 



where 



tq = VLq + glujR, pg 



ilq — nigR 



(4) 



\ng + m*R ' 
iVo-2 = 2R'f,{x)[ng{ng - 1) + mlR/2]/x' , (5) 



with fig = Jx'^ + {m*Ry, m* = niq — g^a, R the bag radius, and Xq the quark Pauli spinor. 
The eigen-frequency, x{xo), of this lowest mode in medium (free space) is determined by 
the boundary condition at the bag surface, 

jo{x) = (3gji{x). (6) 
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The form of the quark wave function in Eq. (j^) is almost identical to that of the solution in 
free space. However the parameters in the expression have been substantially modified by 
the surrounding nuclear medium. Note that as the value of g^W is usually larger than m^, 
the quantity, /3q, becomes larger than unity, which means the lower component of the Dirac 
spinor is enhanced. Thus the quarks in the nucleon embedded in the nuclear medium are 
more relativistic than those in a free nucleon. 

The mean values of the scalar {a) and vector iw) fields in symmetric nuclear matter are 
self-consistently determined by solving the following set of equations: 

(7) 

''d^k^^^ (8) 

(9) 

dR - (1°) 

where p {ps) is the baryon (scalar) density and kp is the nucleon Fermi momentum, g„ = 
3g^S{0), g^j = Sg^, and the quantity, C(w), is defined by 

C{a) ^ S{a)/S{0) = - (^^) /^., (H) 

with S{a) = J\ys,gd^rq{f)q{r). Using the quark wave function for the MIT bag, Eq.(^, S(a) 
can be explicitly evaluated: 

S{a) = [n/2 + mlR{n - l)]/[fi(fi - 1) + mlR/2]. (12) 

The second term in Eq.(P), {—zq/R), parametrizes the sum of the center-of-mass (cm.) 
motion and gluon corrections [|13|. It will be important to note that this center-of-mass 
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treatment is different from that of Jin and Jennings |T^ when we consider later the reduction 
of the bag constant, B. 

For a free nucleon, the quantities, zq and Bq, are determined simultaneously by requiring 
the free nucleon mass to be m,N{p = 0) = 939 MeV and by the stability condition, Eq.(|I^), 
for a given bag radius, Rq, in free space (treated as an input parameter). 



It is worthwhile to note that the self-consistency condition in QMC is identical to that in 
QHD, except that in QHD one has C(w) = 10, which corresponds to a point-like nucleon. 
All information on the internal structure of the nucleon is contained in C{a). 

Equivalently, the equation for the mean field, a, can also be derived by the thermody- 
namic condition ^^^ = 0, where the total energy per nucleon, Etota.\/A, for symmetric 
nuclear matter is given by 

^*°*''- ^ ^'"d'kJm%'{a) + P + ^a^ + -^p. (13) 



A {2'nYpJ V Jv V / i i 2p 2ml 

In the following, we will consider the modified QMC model which was proposed by Jin 
and Jennings [Q, and includes the effect of the reduction in the bag constant, B, (instead 
of a fixed value, Bq) in the nuclear medium. 

The main difference of the approach used in the present work from that of Jin and Jen- 
nings |l^ (see also Fleck et al. 0] and our earlier work) is that we do not compensate for the 
spurious center-of-mass motion by subtracting the average value of the momentum squared 
of the three quarks. It has been shown that this subtraction, which gives a strong field 
dependence, would cause significant reduction of the vector potential. Actual calculations 
using a relativistic oscillator in an external field, where the spurious center-of-mass motion 
can be removed exactly, show that this field dependence is quite small and thus the simple 



subtraction procedure is now believed to be inappropriate [^ . 

When a nucleon is put into the nuclear medium, the bag as a whole reacts to the new 
environment, thus both zq and B may be modified somehow. There has been some con- 



sideration of how Zq might vary depending on the type of baryons p2[ and on the density 



T0| . However, in the present work, we assume that Zq is a fixed constant. In the MIT bag, 
the bag constant, i?o, is introduced to satisfy energy-momentum conservation and provides 
the necessary pressure to confine the quarks. It enters the nonlinear boundary condition, 
— \n^d^^{x)q{x)] = Bq, of the model, where n^ is a unit four vector normal to the bag 
surface. In free space, the bag constant term contributes approximately 200 — 300 MeV to 
the nucleon energy. However, it is expected to be reduced at finite density and temperature. 
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At a sufficiently high density, once nucleons overlap severely, the whole concept of the bag 
is probably no longer realistic. One would expect that deconfinement of quarks could occur, 
and the bag pressure might eventually become zero |2^ . Then one can expect that a decrease 
of the bag pressure should lead to a larger bag radius and hence larger r.m.s. electromagnetic 
radii. On the other hand, these have been constrained by experiment |jT^|-|I^. 

In our previous work [^] we kept zq and Bq fixed and used Eq. ([10|) to determine the 



in-medium bag radius, R, for a finite nuclear density, p. Typically the quark r.m.s. radius, 
r*, calculated by the bag wave function is slightly increased, although the bag radius, R, 
typically decreases by a few percent at normal nuclear density. The extent of nucleon swelling 
in QMC, with a constant value for the bag pressure, Bq, has been shown to be consistent 



with the experimental constraints [^ 



Two models for the reduction in the bag constant, motivated by the work of Jin and 



Jennings |14|, are examined in the following. In the direct coupling model, the bag constant 



is related to the mean scalar field, through 

-Do V niN / 

where g^ and 6 are real, positive parameters. The case, 6 = 1, corresponds to the model 
studied by Blunden and Miller [|l^. In the limit, 6 -^ oo, it becomes a simple exponential 
form, B/Bq = e"^^'' '^Z™'^ . On the other hand, in the scaling model, the bag constant is 
connected to the effective nucleon mass in the medium, 



B frrij^ 



^ \ K. 



(15) 



Bq \mN. 

where k is a real, positive parameter. Note that the standard treatment of the QMC model 
can be recovered in the limits, g^ = 0, for the direct coupling model and, k = 0, for the 
scaling model. 

Because of the additional a dependence arising in Eq.(||), through B(a), the self- 
consistency condition for the scalar field has to be modified. By a straightforward evaluation, 
we may replace the quantity, C(a), in Eq.(|) by, 
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C{a)^Cnc{^)^C{a) + ^g^nR'^(l-^'-^] , (16) 

for the direct coupling model, and 

C{a)^Csc{a)^C{a)(l-^nR'^] , (17) 

for the scaling model. Clearly both models give rise to a smaller bag constant. The values for 
both quantities, Ccci^) ^^^ Cscio'), decrease as the density (the value of the a) increases, 
which is a similar feature to that of the standard QMC model. 

III. ELECTROMAGNETIC RADIUS OF THE NUCLEON 

The evidence concerning a possible nucleon size change in the nuclear medium comes 
largely from lepton-nucleus scattering, in particular, from quasielastic electron-nucleus scat- 
tering. At energy and momentum transfers of several hundred MeV, the virtual photon 
is able to "see" single nucleons. Analyses of the spectrum for the (e, e'N) process at the 
quasielastic peak allow one to extract information for a bound nucleon. The most outstand- 
ing phenomena are perhaps the missing strength in the Coulomb sum rule |T^ and the 



remarkable ?/-scaling behaviour |T^. However, different models often lead to quite different 



interpretations. Calculations with a semiclassical RPA theory suggested that the nucleon 
form factors had to be modified, with the nucleon r.m.s. radii inceased by about 13% in 
^^C, 23% in ^'^Ca and 21% in ^^Fe p5i. This swelling of the proton may be interpreted as 



a distortion of the pionic cloud p6[ or as a reduction of the effective masses of mesons and 



nucleons p7[] . 



The most convincing constraint on the increase of the nucleon r.m.s. radius is from 
the systematic analysis of experimental data covering different momentum transfer regions, 
where an appealing scaling behaviour of the cross section clearly exists. Analogous to 
the Bjorken scaling in electron-nucleon deep-inelastic scattering, this so-called y-scaling in 
quasielastic electron-nucleus scattering reveals the nucleon momentum distribution in the 
nucleus. The scaling feature makes this reaction an ideal testing ground for medium effects 
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on the nucleoli form factors. For recent comprehensive reviews on this subject we refer to 
Ref. [0. 

The idea of y-scahng is based on an analysis under the following basic assumptions ||28|| : 



i) impulse approximation (single nucleon knockout), ii) no final state interactions, and iii) 
nonrelativistic nucleon kinematics. These assumptions can be more or less justified by the 
careful selection of a certain kinematic region. At energies between the giant resonances and 
pion production threshold for intermediate three momentum transfers (e.g., larger than 500 
MeV ) inclusive electron-nucleus scattering is dominated by the quasielastic peak, which 
results from scattering incoherently from individual nucleons embedded in the medium. In 
the laboratory frame, the differential cross section for this process in the plane-wave impluse 



approximation can be expressed as [^,pS 
(Pa 



dVtfdv 



f2|d^k|dEae^iq,uJ■,k,E)S^ik,E) 



*2 



5{u + MA-\jml + k' -^Ml_, + P), (18) 

— * — * — * 

where k and k' = k + q are the initial and final momenta of the struck nucleon, E is the 
energy of the initial nucleon, and (z/, q) is the four-momentum of the virtual photon. S{k, E) 
is the spectral function (the probability of finding a nucleon of momentum k and energy 
E in the nucleus). o"ei is the elementary electron scattering cross section from an off-shell 
nucleon. Ma and Ma-i are the masses of the target and residual nucleus. 
For sufficiently large q, the cross section can be rewritten as 



rf2 



— -ae{q,y)-E{y), (19) 



dfledu 
where q = \q\ and the scaling variable, 

y= (Ma + iy)^jA'^ - Ml_^W^ - gA /W\ (20) 

represents the momentum component of the nucleon parallel to q before the reaction. W = 



[Ma + i^y — (p' is the center-of-mass energy and A = {M\_y — vn\^ + l^^)/2. The scaling 
function, F[y\ is the probability to find a nucleon with momentum fraction y and is a 
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universal function for different combinations of q and v. The quantity, ae{q,y), relates to 
the elementary electron nucleon cross section. 

Note that y-scaling is only valid within a well defined kinematic range. The sensitivity 
of scaling behaviour to the nucleon size can provide a reliable constraint on the possible 
nucleon size change in the medium. Any residual q^ dependence could in principle be 
attributed to modified electromagnetic form factors of the nucleon. In this work, we do not 
try to reproduce the quasi-elastic electron scattering data, instead, we use the best estimate 
of the possible change of size of a bound nucleon available from ^/-scaling data to restrict 
the important model constant — the bag constant, B. 

Although little is known about the off-shell 'yNN vertex, the prescription of de Forest 
||T9| is believed to be successful and is valid at about 5% near the quasielastic peak. By 
ascribing the off-shell correction to a kinematic factor, the single nucleon response is taken 



to be on- mass shell. More general discussions of this subject can be found in Ref. pO[. In 
the rest of this section, we calculate the electromagnetic form factors for an on-shell nucleon 
with the parameters determined in the previous section in a self-consistent fashion. To make 



the predictions as reliable as possible, we use the technique developed in Ref. |]2T]], where 
corrections for spurious center-of-mass motion and relativistic kinematics were considered. 
As is well known, the electric (Ge) and magnetic (Gm) form factors for an on-shell 
nucleon can be most conveniently evaluated in the Breit frame, 

{NAl)\nO)\N,{-l)) =xIxsGe{Q% (21) 

{NAl)\AO)\N,{-p=x'/^XsGMm, (22) 

where Q^ = —q^ = g*^, Xs and xl' ^^^ Pauli spinors for the initial and final nucleons 
respectively. The major advantage of the Breit frame is that Ge and Gm are explicitly 
decoupled and can be determined separately by the time and space components of the 
electromagnetic current (J^)- Note that, in the above definitions, Eqs. ( plf ) and (p2D, both 
the initial and final states are physical states which incorporate meson clouds. 
The electromagnetic current of the quark is simply 
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ji;{x)=Y^Qfeqf{x)rqf{x), (23) 

/ 

where qf{x) is the quark field operator for the fiavor / and Qf is its charge in units of proton 
charge e. The momentum eigenstate of a baryon is constructed by the Peierls-Thouless (PT) 



projection method |21.31 



^PT(fl, X2, X3]p) = A^PTe*^'^"'"g(fi - Xc.m.)q{x2 - Xc.ra.)q{Xz - fern.), (24) 

where A^pt is a normahzation constant, p the total momentum of the baryon, and Xc.m. = 
(xi + X2 + x^)/?> is the center of mass of the baryon (we assume equal mass quarks here). It 
can be shown that the PT wave function satisfies the condition of translational invariance. 
Using Eqs. (|23| ) and (p^, the nucleon electromagnetic form factors for the proton's quark 
core can be easily calculated by 

Ge{Q^) = Jd\UQr)p,{r)K{r)/DpT, (25) 

Gm{Q') = {2m%/Q) Jd'rji{Qr)(3gjo{xr/R)j,{xr/R)K{r)/DpT, (26) 

DpT = ld?rp,{r)K{r), (27) 



where Dpt is the normalization factor, pq{r) = Jq^xt/R) + Pgjf{xr/R), and K(r) = 
j d?u Pqiu) Pq{—u — r) is the recoil function to account for the correlation of the two spectator 
quarks. 

At this stage, there is no satisfactory covariant treatment for the MIT bag model. On the 
other hand, relativistic effects are important for most dynamic variables, especially for form 
factors at large momentum transfer. They lead to a sizable correction for the r.m.s. radius 
of the nucleon. In this paper, we use a semi-phenomenological method to account for the 
relativistic corrections which is consistent with a mean field treatment of the QMC model. 
Since a static MIT bag is an extended spherical object, it would be deformed if it were 
viewed in a moving frame of reference. It is crucial to include this Lorentz contraction of 



the bag for calculating form factors at moderate momentum transfer |]2l| , |32[| . In the prefered 
Breit frame, the photon-quark interaction can be reasonably treated as an instantaneous 
interaction. The resulting form factors can be expressed through a simple rescaling, i.e., 
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GEiQ') = {^?Gt{Q'm%'/E*% (28) 



GMiQ') = i^rGtiQ'm%'/E*% (29) 



where E* = \/m*j^'^ + Q'^/A and G''l^^{Q'^) are the form factors calculated with the spherical 
bag wave function [Eqs. ( p5D and (^)]. The scaling factor in the argument arises from the 
coordinate transformation of the struck quark whereas the prefactor, {m*j^ / E*y , comes from 



the reduction of the integral measure of two spectator quarks in the Breit frame [Q . This 
prefactor is not unambiguously determined, but it is far less important in practice than the 
Lorentz contraction which rescales Q^. The corresponding electromagnetic radii squared are 
thus given by, 



2 ^ / 2\ _ _« dGE,M{Q ) 
rE,M - [r ;^ .^^ - b ^^2 



(30) 



IV. PION CLOUDS AND MESON EXCHANGE CURRENTS 

Before we proceed further it is useful to review the role of meson clouds. As is well- 
known, a realistic picture of the free nucleon should include the surrounding meson clouds. 
Following the cloudy bag model (CBM) PB| , |5^ , we limit our considerations to the meson 



cloud corrections from the most important component, namely the pion cloud. In free 
space, the pion is a Goldstone boson in the chiral limit and is necessary to maintain chiral 
symmetry. In an infinite nuclear medium, the mean value of the pion field vanishes due to 
parity conservation, however, the virtual process of pion exchange is always possible. 

In the QMC model, the medium effect on the nucleon is induced by the a and u meson 
fields. The solution for a bound nucleon has already incorporated certain many body corre- 
lations due to the self-consistent solution of the equations of motion. Estimates show that 
the contributions of meson exchange currents (MEG) to the quasielastic cross section do not 
show the y-scaling feature, and are thus expected to be irrelevant for the present study ^5 . 



One of the major implications of the ^/-scaling analysis of quasielastic scattering data is that 
MEG corrections appear to contribute primarily to the behaviour of the nucleon momentum 
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distribution, rather than to a change of the effective electromagnetic form factors of the 
nucleon. Hence we consider it is reasonable to neglect possible explicit two- or more-body 
MEC contributions for the present study. Thus, we limit the meson cloud corrections to 
one-pion loop on the same nucleon in the usual spirit of the CBM. 

The Lagrangian related to the pion field and its interaction (we use quark-pion pseu- 
doscalar coupling), up to leading order of I/Ztt, is pSJ 



1 1 i 

^^q = ^id^-Kf - i^rnliz'^ - Yfll^r ■ TzqSs, (31) 

where 6s is a surface delta function of the bag, m.,, the pion mass and /tt the pion decay 
constant. The electromagnetic current of the pion is 

j^(x) = -ie['K\x)d^7T{x) - TT{x)df'7r\x)], (32) 

where 7i{x) = -7s[7ri(a;) + i'JT2{x)] either destroys a negatively charged pion or creates a 
positively charged one. As long as the bag radius is above 0.7 fm, the pion field is relatively 
weak and can be treated perturbatively [^]. A physical baryon state, \A), can be expressed 
asg 

\A) ^ ^[1 + {niA -Ho- AHjAy'Hj] \Ao), (33) 

where A is a projection operator which projects out all the components of \A) with at least 
one pion, Hj is the interaction Hamiltonian which describes the process of emission and 
absorption of pions. The matrix elements of Hj between the bare baryon states and their 
properties are [^| 



<Ck) - {Ao\H,\n,ik)Do) = ^ p^f^lp E^S^.^(^:. ■ k)C^TT.tii: ■ e,), (34) 



.DA, 



AD(T\ _ „.AD/ 



<^(fc) = {AoiT,{k)\Hj\Do) = v^^ik) = -v^fik) = v^n-kl (35) 



where the pion has momentum k and isospin projection j, /g is the reduced matrix element 



for the vrDo -^ Aq transition vertex, u{kR) = 3ji{kR)/kR, uok = Jk'^ + m^, and Sm and in 
are spherical unit vectors for spin and isospin, respectively. Note that the nNN form factor 
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u{kR) is fully determined by the model itself and depends only on the bag radius. The 
bare baryon probability in the physical baryon state, Z^, determined by the normalization 
condition, is 



Z^ 



i + Y^llil\' ^^r dkk%\kR) 



^ y m^ / 127r2 Jq uJk{mA - mo - uJkf 
where P denotes a principal value integral. 



(36) 



For the detailed expressions for the pion loop contributions, we refer to Ref. ||21I|. The 
only modifications one should consider, for the purposes of studying nuclear matter, are 
to replace the bare masses and coupling constant with those in medium, i.e., to make the 
substitutions: m^ -^ m* , m^ -^ m*jj and f^AD -^ Had- ■'■^ principle, the existence of the tt 
and A inside the nuclear medium would change the mean fields and probably lead to some 
modification of their properties. Since the pion is well approximated as a Goldstone boson, 
the explicit chiral symmetry breaking is small in free space, and it should be somewhat 
smaller in nuclear medium 0. While the effective pion mass, m*, may be slightly changed 
in the medium |^B|], we still use m*^ = 171^^, because the pion field contributes relatively little 
to the form factors (other than the neutron charge form factor, Geii)- As the A is treated on 
the same footing as the nucleon in the CBM, its mass should vary in a similar manner to that 
of the nucleon. Thus we assume that the in-medium and free space N — A mass splitting are 
approximately equal, i.e., m\ — m*^ ~ ttia — itin- This approximation is exactly satisfied in 
the standard QMC model (See Ref. 0]). In free space, the physical it AD coupling constant is 
obtained by f^^ ~ ( ^S^w ) f^^. There are corrections to the bare coupling constant, f^^ , 
such as those from the nonzero quark mass and the correction from spurious center of mass 
motion. Therefore, following the standard treatment, in practice we use the renormalized 
coupling constant in our calculation here, /^^ ^ 3.03, which corresponds to the usual ttNN 
coupling constant, /^^vat — 0.081. In the medium, the vrA^A^ coupling constant might be 
expected to decrease slightly due to the enhancement of the lower component of the quark 
wave function, but we shall ignore this small density dependence in the present study and 
use the approximate relation, f*j^j^ — f^NN, for simplicity. 
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V. DISCUSSION OF THE RESULTS 

In the previous standard (i.e., "constant -Bq") QMC models, typical solutions show that 
the values of the scalar mean field are quite large - the term g^W is of the order of 100 MeV 
at normal nuclear density, and thus the final results are insensitive to the values used for the 
current quark mass, niq P,p!3| . Later variations introduced some density dependence of the 
bag constant, B, through a parametrized, phenomenological form (static parameter g^ and 5 



or k) \TA\. Thus for a given bag radius (input for free nucleon), there are two truly dynamical 
free parameters, the quark and meson coupling constants, g^ and g^^, which are adjusted 
to reproduce the saturation properties of symmetric nuclear matter. The effective nucleon 
mass, the bag radius in matter, and the value for the scalar mean field are all determined 
by the equations of motion, whereas the mean field of the vector meson is directly related 
to the baryon density (see also Eqs.(0) - (|l3)). A typical set of parameters and the value of 
the nuclear incompressibility, K, obtained for each model - direct coupling models I and II 
(DC-I and DC-II), the scaling model (SCALE) and the standard QMC - are listed in Table |. 
The properties of these self-consistent solutions as a function of the nuclear density (in units 
of normal nuclear matter density, po = 0.15 fm^'^), are shown in Figs. 0-|^, where the bag 
radius in free space is taken to be 0.8 fm. For a study of the (rather weak) dependence of 
the final results on the values chosen for the bag radius in free space, see Refs. P,p!3|]. 

In Fig. ID, r* (r^) denotes the quark r.m.s. radius calculated using the quark wave function 
in medium (free space), and m^ {ttin) denotes the nucleon mass in medium (free space). The 
density dependence of the effective nucleon mass is more or less similar for all the QMC type 
models considered in this study, and the reduction of the mass is very moderate compared 
with that of QHD. However, the QMC type models among themselves show a quite different 
density dependence of the quark r.m.s. radius, r*. 

In Fig. 1^, we show the energy per nucleon (with respect to the free nucleon mass, E/A — 
771 at) calculated for each model, together with the the value of the nuclear incompressibility, 
K. The parameters in each model are fitted to reproduce the energy at the saturation point, 



17 



—15.7 MeV. All the QMC type models reproduce reasonable values (200 ~ 300 MeV) for 
the nuclear incompressibility, K (for the parameter sets listed in Table |I|). As is well known, 
the value calculated in a naive version of the QHD model is too large, and thus, the energy 
per nucleon grows more rapidly than those of the QMC type models as the density increases. 

We also show the density dependence of the bag constant, -B, and bag radius, i?, in 
Fig. ||. Although the bag radius calculated for the standard QMC model decreases slightly 
as the density increases, the increase of the bag radius calculated for the modified QMC 
models is very pronounced. As expected, a large reduction of the bag constant causes a 
large increase of the bag radius. 

The quality of the nucleon electromagnetic form factors in free space has been described 
in our previous publication [^]. The corrections for the center-of-mass motion and Lorentz 



contraction lead to a significant improvement over earlier static bag model calculations, in 
particular, at moderate momentum transfers. The pionic correction tends to soften the both 
form factors, Ge and Gm-, and thus increase the corresponding electromagnetic r.m.s. radii 
of the proton, Te = (r\) and Tm = (^m) • (The corresponding quantities in free space 
will be denoted by r^ and r%, respectively.) We also note that the size of these corrections 
is larger when smaller values of the bag radius in free space are chosen. 

As in our previous study [|^, the charge form factors are much more sensitive to the 



nuclear medium than the magnetic ones. The variation of the magnetic form factors with 
respect to their free space counterparts is nearly one order of magnitude smaller than that for 
the charge form factors. As for the density dependence of the electromagnetic form factors, 
they are clearly suppressed at small Q^, relative to those in free space. This indicates 
that the electromagnetic radii of the bound nucleon have increased. In other words, the 
nucleon becomes "swollen" once it is put into the nuclear medium. For a fixed Q^ (less than 
0.3 GeV^), the form factors decrease almost linearly with the increase in the nuclear density. 
At Q^ ~ 0.3 GeV^, the proton and neutron charge form factors are reduced by roughly 5% 
and 6% at p = 0.5po, and 8% at normal nuclear density, po = 0.15 fm^^. Similarly, the 
proton and neutron magnetic form factors are 1% and 0.6% smaller at p = 0.5po, and 1.5% 



and 0.9% at normal nuclear density. 

As explained earlier, the best experimental constraints on the changes in these electro- 
magnetic form factors come from the analysis of y-scaling data. In the kinematic range 
normally covered by this analysis, the electron-nucleon scattering cross section is predom- 
inantly magnetic, so that the limit extracted applies essentially to Gm- For example, in 
Fe it was concluded by Sick that the nucleon r.m.s. radius cannot vary by more than 3% 
T6| . For the standard QMC model, with a fixed bag constant, the calculated increase in the 



r.m.s. radius of the magnetic form factors is less than 0.8% at po- This is well within the 
experimental limit. For the modified QMC models (with smaller B), however, the variation 
of the r.m.s. radius can be much stronger. 

In Tables || — 0, we show the changes of the electromagnetic r.m.s. radii of the proton, 
Ste^m = ijE,M — f% m)/''^e Mi ^^ ^^c bag constant for three different nuclear densities. Two 
values of the bag radius in free space are used in this calculation. In the direct coupling 
model, we take the limit, 5 ^ oo, and vary g^ . For each selected model (fixed g^ or 
k), the coupling constants, gl/^ii and 5'^/47r, are first determined by fitting the saturation 
properties of nuclear matter at normal nuclear density. The mean fields, quark frequency, 
and effective nucleon mass are then subsequently solved self-consistently for other nuclear 
densities. The current quark mass is taken to be zero. It is clear from the Tables that the 
bag constant and the electromagnetic r.m.s. radii are strongly correlated. A smaller bag 
constant will lead to a larger size of the bound nucleon. 

For finite nuclei, it is well known that a large portion of the nucleons are distributed 
near the nuclear surface. The experimental constraints thus correspond to the quantities at 
the average nuclear density, which we take to be about O.TpoQ- Fig- @ shows the correlation 
between the bag constant and the magnetic radius of the proton at a nuclear density p = 



^Using standard nuclear density distributions the average densities of ^^O, '^^Ca and ^^Fe are 
0.6po; 0.7/90 and 0.7 po, respectively. 
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0.7po- For an increase of about 3% in the magnetic radius {Svm — 3%), which is the upper 
hmit allowed from the y-scaling data, the bag constant is allowed to decrease by about 8% 
for the direct coupling model and 16% for the scaling model, with the bare bag radius, 
Rq = 1.0 fm. These reductions become 12% and 18% with a smaller bag radius, Rq = 0.8 
fm. 

The possible reduction of the bag constant, B, might also be extracted from the electric 
form factors if we had a reliable experimental constraint for the charge radius. In the y- 
scaling analysis, since the electric and magnetic form factors contribute typically in the ratio 
1:3 in the kinematic regime selected, the corresponding limit on Ge would be nearly 10%. 
For the electric form factors the best experimental limit seems to come from the Coulomb 



sum-rule, where a variation bigger than 4% would be excluded ||T8[. With this strong limit, 
the bag constant is allowed to be reduced by roughly 10% in the scaling model and is almost 
unchangeable in the direct-coupling model . 

VI. CONCLUSION 

In summary, the bag constant in the nuclear medium, as a crude representation of the 
very complicated QCD mechanism of confinement, may be expected to decrease. We have 
examined the possible variations of the bag constant in two modified versions of the QMC 
model. The parameters are determined to reproduce the empirical values for the bulk 
properties of symmetric nuclear matter. We have calculated the density dependence of 
the electromagnetic form factors of the bound nucleon in the medium, then compared the 
corresponding r.m.s. radii with the best available y-scaling analysis for quasielastic electron- 
nucleus scattering. Since the change in the size of the nucleon is strongly limited by these 
y-scaling data, so is the bag constant. Our results show that the bag constant in medium, 
B, is allowed to decline by only 10-17% at the typical average nuclear density of about 
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0.7poQ- This limit would be reduced by a further 7% at normal nuclear density. This is a 



considerably smaller change than most of the cases studied by Jin and Jennings [Q, but 
comparable with the usual estimates of the change in the gluon condensate [p^ , or of the 



bag constant itself |^ in nuclear matter. 

It is worth mentioning that the pion cloud of the nucleon plays an important role in 
calculating the electromagnetic form factors. Without the pion cloud, the nucleon r.m.s. 
radius tends to be somewhat small, for a reasonable bag radius. Thus the relative variation 
is much more sensitive to the bag constant. Typically, 5rM without the pionic correction 
is three times larger than that with the pion field. We would like to emphasize that the 
MEC effects which can be important in conventional nuclear physics have not been included 
in the present study. This should be reasonable in a first treatment because MEC correc- 
tions are not very important near the quasielastic peak and they contribute primarily to the 
behaviour of the nucleon momentum distribution, rather than to a change of the effective 
electromagnetic form factors of the nucleon. Further improvements rely on a better treat- 
ment of the QMC model and require the explicit inclusion of nucleon-nucleon correlations 
and medium polarization effects. 

We are pleased to thank I. Sick for helpful correspondence concerning the constraints 
from y-scaling analysis. We are also grateful to M. Ericson for informative discussions. This 
work was supported by the Australian Research Council. A.W.T and K.S. acknowledge 
support from the Japan Society for the Promotion of Science. 



^In terms of the actual models studied this means that g^ should be less than 1.2 in the direct 
coupling model with 5 ^ oo (c.f. Eq. (|l^)) and k less than 1.2 in the scaling model (c.f. Eq. (|l5|)). 
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TABLES 
TABLE L Saturation properties calculated in each model for symmetric nuclear matter at 

normal nuclear density, po = 0.15 fm^^. DC-I and DC-II refer to the direct coupling models with 

g^ = A and 5 = 890, and g^ = 2.8 and 6 -^ oo, respectively. SCALE denotes the scaling model 

with K = 1.45 and QMC denotes the standard QMC model treatment. (The bag radius in free 

space is taken to be Rq = 0.8 fm.) x (xq) is the bag eigenfrequency in matter (vacuum), and K is 

the nuclear incompressibility. 



model 


5^/4^ 


5^/4^ 


x/xq 


m^/niN 


R/Ro 


B/Bo 


K{MeV) 


DCT 


0.521 


1.600 


0.9393 


0.8380 


1.239 


0.636 


278.2 


DC-II 


2.811 


6.393 


0.8354 


0.7557 


1.160 


0.707 


242.3 


SCALE 


3.185 


8.525 


0.8941 


0.7907 


1.079 


0.711 


339.0 


QMC 


5.373 


5.257 


0.8368 


0.8048 


0.982 


1.000 


278.2 



TABLE II. Variation of the bag constant vs the electromagnetic radii of the proton in the 
direct coupling model with the limit, 6 -^ oo, for Rq = 0.8 fm. The calculated values of the 
electromagnetic r.m.s. radii in free space are, r^ = {r'^)^ = 0.736 fm and r^ = {r"^) ^j = 0.731 
fm. Below, 5rE,M = {fE,M — r% m)/''"e m ^^'^ ^^^ given as percentages. 









p/f 


= 0.5 




p/p 


= 0.7 




p/p 


= 1.0 




9^ 


5^/4^ 


ffS/4vr 


B/Bo 


6rE 


Stm 


B/Bo 


6rE 


^rM 


B/Bo 


6rE 


SrM 


0.4 


5.015 


5.604 


0.976 


3.20 


0.82 


0.968 


4.26 


1.10 


0.958 


5.66 


1.51 


0.8 


4.651 


5.906 


0.952 


3.86 


1.37 


0.937 


5.17 


1.85 


0.916 


6.93 


2.52 


1.2 


4.293 


6.168 


0.928 


4.56 


1.95 


0.905 


6.16 


2.64 


0.874 


8.32 


3.61 


1.6 


3.935 


6.373 


0.903 


5.30 


2.56 


0.872 


7.22 


3.48 


0.832 


9.84 


4.78 


2.0 


3.575 


6.502 


0.879 


6.08 


3.19 


0.840 


8.35 


4.37 


0.790 


11.5 


6.03 
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TABLE III. Same as in Table || except Rq = 1.0 fm. The free space values are now, r^ = 0.83 



fm and r^^ = 0.80 fm. 









p/p 


= 0.5 




p/p 


= 0.7 




p/p 


= 1.0 




9^ 


5a/4^ 


ffS/4vr 


B/Bo 


SrE 


SrM 


B/Bo 


6rE 


SrM 


B/Bo 


6rE 


6rM 


0.4 


4.692 


4.837 


0.977 


3.81 


1.07 


0.970 


4.99 


1.41 


0.960 


5.78 


1.16 


0.8 


4.341 


5.140 


0.954 


4.60 


1.65 


0.940 


6.07 


2.20 


0.921 


7.95 


2.92 


1.2 


4.001 


5.418 


0.931 


5.45 


2.28 


0.909 


7.24 


3.05 


0.880 


9.58 


4.10 


1.6 


3.668 


5.657 


0.907 


6.35 


2.95 


0.878 


8.51 


3.98 


0.840 


11.4 


5.39 


2.0 


3.340 


5.844 


0.883 


7.31 


3.66 


0.846 


9.89 


4.99 


0.799 


13.3 


6.81 



TABLE IV. Variation of the bag constant vs the electromagnetic radii of the proton in the 
scaling model with Rq = 0.8 fm. The values of the electromagnetic r.m.s. radii in free space are, 
r^ = 0.74 fm and r^j = 0.73 fm. See also the caption of Table O. 









p/p 


= 0.5 




p/p 


= 0.7 




p/p 


= 1.0 




K 


5^/471" 


95/4^ 


B/Bo 


6rE 


SrM 


B/Bo 


6rE 


6rM 


B/Bo 


6rE 


SVM 


1.0 


3.732 


7.219 


0.888 


3.75 


1.73 


0.850 


4.70 


2.11 


0.800 


7.21 


3.43 


1.2 


3.468 


7.727 


0.866 


4.02 


2.03 


0.821 


5.59 


2.84 


0.761 


7.88 


4.07 


1.4 


3.241 


8.359 


0.843 


4.34 


2.35 


0.791 


6.08 


3.32 


0.721 


8.67 


4.79 


1.6 


3.021 


9.058 


0.819 


4.67 


2.68 


0.759 


6.63 


3.83 


0.681 


9.56 


5.57 


1.8 


2.652 


9.215 


0.805 


4.78 


2.89 


0.740 


6.81 


4.13 


0.654 


9.95 


6.06 



26 



TABLE V. Same as in Table |^ except Rq = 1.0 fm. The values of the electromagnetic r.m.s. 
radii in free space are, r^ = 0.83 fm and r^j = 0.80 fm. 









P/P 


= 0.5 




p/p 


= 0.7 




p/p 


= 1-0 




K 


9^/4vr 


5S/4vr 


B/Bo 


6rE 


SVM 


B/Bo 


6rE 


SVM 


B/Bo 


6rE 


(Jr-Af 


1.0 


3.431 


6.125 


0.899 


4.27 


1.87 


0.865 


5.80 


2.58 


0.822 


7.91 


3.59 


1.2 


3.183 


6.565 


0.879 


4.56 


2.16 


0.839 


6.23 


2.99 


0.787 


8.58 


4.21 


1.4 


2.952 


7.064 


0.859 


4.87 


2.46 


0.813 


6.71 


3.44 


0.752 


9.33 


4.88 


1.6 


2.791 


7.814 


0.835 


5.31 


2.84 


0.781 


7.36 


3.98 


0.712 


10.38 


5.73 


1.67 


2.639 


7.762 


0.832 


5.30 


2.87 


0.777 


7.35 


4.04 


0.705 


10.40 


5.83 
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FIG. 1. Density dependence of the effective nucleon mass and quark r.ni.s. radius. The param- 
eters used in the calculation are given in Table |. 
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FIG. 2. Calculated values of the energy per nucleon (with the free nucleon mass subtracted) 
for symmetric nuclear matter. The values, K^ denote the nuclear incompressibility calculated in 
each model. 
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FIG. 3. Density dependence of the bag constant and the bag radius. 
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FIG. 4. Correlation of the bag constant with the magnetic radius of the proton at the average 
nuclear density, p/po = 0.7. The number in the bracket refers to the bare bag radius, Rq. 
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